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We study ground state properties of periodic Anderson model in a two-dimensional square 
lattice with variational Monte Carlo method. It is shown that there are two different types 
of quantum phase transition: a conventional antiferromagnetic transition and a Fermi-surface 
reconstruction which accompanies a change of topology of the Fermi surface. The former is 
induced by a simple back-folding of the Fermi surface while the latter is induced by localization 
of / electrons. The mechanism of these transitions and the relation to the recent experiments 
on Fermi surface are discussed in detail. 
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1. Introduction 

Recently, antiferromagnetic (AF) transition which ac- 
companies an abrupt change of Fermi surface has been 
widely observed in heavy fermion compounds and has at- 
tracted much attention. De Haas- van Alphen (dHvA) ex- 
periments in Ce compounds, such as CeRh2Si2,^ Celns^ 
and CeRhlns,'^ show that the topology of Fermi surface 
abruptly changes across the AF transition under pres- 
sure. According to band calculations, Fermi surface in 
the AF phase is well explained by a 4/-localized model 
while that in the paramagnetic phase is well explained 
by a 4/-itinerant model. This Fermi-surface change 
has also been suggested in the Hall effect experiment in 
YbRh2Si2.^ Although the Hall coefficient shows a contin- 
uous change as a function of magnetic field at finite tem- 
perature, its slope becomes sharper with decreasing tem- 
perature. From this, it is claimed that there will be a dis- 
continuous jump of the Hall coefficient at the AF quan- 
tum critical point (QCP) at T = 0. This means that a 
change from a "small Fermi surface" (AF side) to a "large 
Fermi surface" (paramagnetic side) occurs at QCP. In 
these experiments, the AF transition and the abrupt 
change of the Fermi surface seem to occur simultane- 
ously. On the other hand, in CeRhi_j;Coa;In5, an abrupt 
change of the Fermi surface is observed at x = 0.3 0.4 
while the AF QCP exists at a; = 0.7 ~ 0.8.*^''^ In this 
case, the AF transition and the abrupt change of the 
Fermi surface are not necessarily coincident and the re- 
lation between them has not been clarified yet. 

In a conventional picture of the QCP in heavy fermion 
system, the AF transition is induced by the competition 
between Kondo screening and RKKY interaction.^ In 
this case, the conduction electrons (c electrons) and the 
localized / electrons hybridize with each other and form 
the composite quasiparticles. The / electrons obtain the 
itinerancy through the hybridization and the AF transi- 
tion is expected to be described by a so-called Moriya- 
Hertz-Millis (MHM) theory.^"" This picture have been 
extensively studied so far and succeeded in describing 
the unconventional phenomena such as a non-Fermi- 
liquid behavior around the QCP. However, it is not clear 
whether the MHM theory can be applied or not to the 



case where the / electrons are strongly localized and 
the itinerant picture becomes inappropriate. Indeed, the 
change of the topology of the Fermi surface mentioned 
above cannot be explained by the simple back-folding of 
the Fermi surface which is based on the itinerant picture. 

Recently, on the other hand, novel type of theory which 
considers the localization of / electrons has been pro- 
posed by some authors. ^^"^^ It is based on the breakdown 
of the Kondo screening or the decoupling of c and / elec- 
trons. In this case, the change of nature of the / electrons 
at the QCP has been emphasized, but the details are still 
uncertain. 

With these facts in mind, we studied the ground state 
phase diagram of the Kondo lattice model (KLM) in 
the previous work.^^ This model is one of the simplest 
models which contain c and / electrons and describe the 
heavy-fermion system. We have found that there are two 
types of quantum phase transitions: a conventional AF 
transition and a topological transition of a Fermi surface 
(Fermi-surface reconstruction). The former transition is 
induced by the folding of the Fermi surface and therefore 
we denote as "conventional" . In the latter transition, on 
the other hand, the Kondo screening becomes almost ir- 
relevant (but does not disappear) and almost fully po- 
larized AF state appears accompanied with the change 
of the topology of the Fermi surface. 

In this paper, we study the ground state properties of 
the periodic Anderson model (PAM). The KLM will not 
be enough to discuss the issue of the itinerancy and lo- 
calization of / electrons, since the / electrons are always 
localized in the KLM. On the other hand, in the PAM, 
we can control the number of / electrons and can discuss 
the change of the nature of / electrons by changing the 
parameters such as U (on-site Coulomb interaction), V 
{c-f hybridization) and Ef (energy level of / electrons). 
Moreover, we discuss the properties of the Fermi sur- 
face in detail by calculating the momentum distribution 
function which was not discussed in the previous paper 
on the KLM.^^ We show that the Fermi-surface recon- 
struction observed in the KLM also occurs in the PAM 
under the condition of U — oo and sufficiently low val- 
ues of Ef. The Fermi-surface reconstruction have two re- 
markable aspects as a quantum phase transition. It can 
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be regarded as a kind of Lifshitz transition, where the 
topology of the Fermi surface changes from hole-Uke to 
electron- Hke. At the same time, it can be regarded as 
a kind of Mott transition, where the / electrons change 
their character from itinerant to localized. 

The organization of this paper is as follows. In §2, we 
introduce the model and variational wave functions. In 
§3, we show the results of VMC calculation. The ground 
state phase diagrams and some physical quantities are 
studied in detail. In §4, we discuss the relation between 
our results and the experiments. §5 is devoted to the 
summary. 

2. Model and Method 

We study the following PAM in a two-dimensional 
square lattice, 



ia i 

where e^, Ef, V and U represent energy dispersion of c 
electrons, energy level of / electrons, hybridization be- 
tween c and / electrons and on-site Coulomb interaction 
between / electrons, respectively. In the PAM, the num- 
ber of / electrons at each site is not fixed and therefore 
the total electron density n = ric + rif is a. controlling 
parameter. 

The trial wave function for a paramagnetic state is 
expressed as, 

I*) = Pi 1$) = Pi nK4. + vkfi) io) (2) 
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Here Uk {vk) corresponds to the mixing amplitude of 
c (/) electrons. Pq is a Gutzwiller projection operator 
which reduces the probability of double occupancy of / 
electrons at the same site. Since we consider the case of 
U = oo, we set g — m the following. The improjected 
wave function |$) is obtained by diagonalizing the fol- 
lowing one-body Hamiltonian. 
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It is just the f7 = PAM with "effective" hybridiza- 
tion V and "effective" /-electron level Ef.V and Ef are 
variational parameters and they are optimized so as to 
minimize the variational energy. This form of trial wave 
function is also used in ref. 18 for a one-dimensional case. 



As an AF trial wave function, we construct a state |$) 
by diagonalizing the following Hamiltonian, 
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with Q = (tt, tt) being an AF ordering vector. Here Ac 
and Af represent the variational parameters of AF gap 
in c and / electrons, respectively. The summation over k 
is carried out in the folded AF Brillouin zone. If we set 
Ac = A/ = 0, eq. (7) reduces to eq. (6). Although the 
choice of the basis and the notation of variational param- 
eters are slightly different from those used in ref. 17, the 
obtained unprojected band dispersion is the same with 
the previous paper^'' (A / and Ac corresponds to M and 
m in ref. 17). 

We optimize the four variational parameters 
{V,Ef,Af,Ac) and find the lowest energy state in 
a V-n phase diagram. Possible candidates for the 
ground state are classified according to the shape of the 
Fermi surface and the band dispersion of the one-body 
part |$): paramagnetic metal (PM), AF metal with hole- 
like Fermi surface (AFh), AF metal with electron-like 
Fermi surface (AFo) and AF metal with V ^ (AFs, S 
denotes "small"). This notation is the same with that in 
the previous paper^'' and the schematic Fermi surfaces 
are shown in Fig. 1. In PM, AFh and AFe, V is finite 
and the one-body part forms the c-f hybridized band. 
It is based on the concept of "large" Fermi surface. In 
contrast, = in AFg, and c and / electrons form 
their band dispersions individually. AFs is based on the 
concept of "small" Fermi surface. However, note that 
the expressions of "large" and "small" Fermi surface are 
not appropriate in the AF state since the volumes of 
the Fermi surface are the same in AFh, AFe and AFs. 
Moreover, AFo and AFg have the same Fermi-surface 
topology and cannot be distinguished from each other 
only by the Fermi surface. The difference between them 
are discussed in the next section. 

3. Results 

3.1 Ground St,ate Phase Diagram 

We show the ground state phase diagrams for different 
values of Ef in Fig. 1. The system sizes are set to be 8x8, 
10x10 and 12x12 with periodic-antipcriodic boundary 
condition. When Ef = 0.0, the state changes from PM 
to AFh Sit V = Vaf as shown in Fig. 1(a). The first 
Brillouin zone is folded and the AF gap opens along 
the newly formed AF Brillouin zone boundary. This is 
a "conventional" second-order AF transition. The region 
of AFh extends to V ^ 0. On the other hand, when 
Ef = —0.2, a novel type of quantum phase transition 
occurs and AFe is realized for V < Vps, as shown in 
Fig 1(b). At V = Vpg, the band dispersion discontinu- 
ously changes from convex upward (AFh side) to convex 
downward (AFe side). It is regarded as a kind of Lif- 
shitz transition, where the topology of the Fermi surface 
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changes. We call this transition a "Fermi-surface recon- 
struction" and its origin is discussed in §3.2 in detail. 
As Ef becomes lower, the region of AFe becomes larger 
and finally overlaps the region of AFh. As a result, when 
Ef = -1.0 (Fig. 1(c)), the direct transition from PM to 
AFe appears for n < 1.78. We can see that the stability 
of AFe greatly depends on the value of Ef. 

For fixed values of Ef and n, the change of the ground 
state is classified into three cases: (i) PM — » AFh, (h) PM 

AFh AFe and (iii) PM AFe. In Fig. 2, we show 
the energy comparison between the trial wave functions 
for each case and discuss the behaviors of the staggered 
magnetization of / electrons, Mf = (rij^.^) — (tIaj.) — 

~{{^Bi) ~ ('^Bj,))' and the average /-electron density 
(n/). 

Figure 2(a) shows the case of (i) at Ef = 0.0 and 
n = 1.880. The groimd state is PM for V > Vaf and AFh 
for V < Vaf- The variational energy of AFh smoothly 
merges with that of PM (horizontal axis) a.t V = Vaf, 
indicating the second-order (continuous) AF transition. 
Mf continuously increases from zero as V decreases from 
Vaf- {n-f) monotonically increases as V decreases from 
the PM side, and shows continuous change across the AF 
transition. 

The case of (ii) at Ef = —0.2 and n = 1.880 is shown 
in Fig. 2(b). As in the case of (i), the transition from PM 
to AFh &t V = Vaf is second-order. The changes of Mf 
and (n/) across the AF transition are continuous. On 
the other hand, when the Fermi-surface reconstruction 
(AFh — AFe) occurs, abrupt changes are observed. At 

V = Vfs, the energy crossing occurs and the ground 
state changes from AFh to AFc. It is clearly a first-order 
transition and physical quantities such as Mf and (n/) 
discontinuously increase across the transition. As one can 
see from Fig. 2(b), the / electrons are almost localized 
in AFe and Mf ~ (n/) is satisfied. It shows that the 
Kondo screening, which works so as to suppress the AF 
order, becomes almost irrelevant in AFc and the local 
/-electron moment is almost fully polarized. 

The case of (iii) at Ef = —1.0 and n = 1.750 is shown 
in Fig. 2(c). Contrary to the cases of (i) and (ii), AFh 
does not appear for any values of V; the variational en- 
ergy of AFh is always higher than that of AFe, and 
the ground state directly changes from PM to AFp at 

V = Vaf+fs- It is a first-order transition, where the AF 
transition and the Fermi-surface reconstruction occur si- 
multaneously. As a result, Mf shows discontinuous jump 
from zero to a finite value a,t V = Vaf+fs- As in the 
case of (ii), Mf (n/) is satisfied and the almost fully 
polarized AF state is realized in AFc. 

Note that AFs does not appear for all cases. It is a 
natural result since AFs cannot gain the mixing energy 
because of V^ = and it is energetically unfavorable un- 
der the condition of y ^ 0. The introduction of small but 
finite value of V does not change the shape of the Fermi 
surface very miich. but it leads to a mixing energy gain 
and stabilizes AFc. In AFc, the c electrons enlarge the 
occupied phase space and become heavy compared with 
the case of AFs bccaiise of the finite value of V . Namely, 
AFe and AFs have almost the same Fermi surface but 



their internal structures of quasiparticles are completely 
different. 

3.2 Origin of the Fermi Surface Reconstruction 

In §3.1, we have shown the existence of the novel type 
of quantum phase transition of "Fermi- surface recon- 
struction" . In this section, we discuss its origin. Figure 3 
shows the band dispersions of the three different states 
obtained by diagonalizing the one-body part |$) with op- 
timized parameters: PM (F=1.50), AFh (F=1.00) and 
AFe (V=0.70) at Ep = -1-0 and n = 1.880. In PM, only 
the lower hybridized band is shown since the upper one 
is away from the Fermi energy and mioccupied. In AFh, 
the energy band is folded and the AF gap opens along 
the AF Brillouin zone boundary ((tt, 0) and (7r/2,7r/2) 
in Fig. 3(b)). This folding docs not change the topology 
of the Fermi surface compared with Fig. 3(a). As V de- 
creases in AFh, the optimized value of V also decreases 
and that of A/ increases. Following the growth of the 
AF gap, the band dispersion changes in the direction 
as shown by the arrows in Fig. 3(b) and finally changes 
into a dispersion which is convex downward. As a result, 
quasiparticle occupation changes in the fe-space and the 
Fermi surface is reconstructed from hole-like (AFh) to 
electron-like (AFe) at V = Vfs- According to this mech- 
anism, it is in principle possible that the Fermi-surface 
reconstruction becomes a second-order transition. How- 
ever, within our trial wave function, the Fermi-surface 
reconstruction is clearly a first-order transition which ac- 
companies the energy crossing. 

Next we discuss the origin of the Fermi-surface re- 
construction from a viewpoint of energy gain. The vari- 
ational parameter V corresponds to the effective hy- 
bridization and thus it will be directly related to the 
degree of Kondo screening. It works so as to gain the 
mixing energy ii^niix, the third term in eq.(l). On the 
other hand, A y corresponds to the AF gap and controls 
the magnitude of the AF order. It stabilizes the local /- 
electron moment and suppresses the hybridization, lead- 
ing to gain the kinetic energy of c electrons i?kin and the 
/-electrons energy En^ , which are the first and the sec- 
ond terms in eq.(l), respectively. In PM and AFh, the 
role of V is dominant and the system forms the hole-like 
Fermi surface which is favorable for i^mix- In AFe, on the 
other hand, the role of A/ exceeds that of V and the 
electron-like Fermi surface, which is favorable for -Ekin 
and Enf, is formed. Namely, the origin of the Fermi- 
surface reconstruction is a change of the mechanism of 
the energy gain. In AFe, the localization of / electrons is 
essential for the energy gain and the AF order is consid- 
ered to be additional. Therefore, the region of AFe in the 
phase diagram mainly depends on Ef, but not so much 
on the total electron density n which is closely related to 
the stability of the AF order. It is in contrast to the case 
of AFh whose stability greatly depends on n. 

The existence of the Fermi-surface reconstruction is 
also suggested in the recent study of the cellular dynam- 
ical mean-field theory (CDMFT) by De Leo et al.'^^ They 
have calculated spectral functions for paramagnetic and 
AF states and shown that the AF transition is accompa- 
nied by a dramatic rearrangement of the spectral weight. 
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Fig. 1 . Ground state phase diagram of (7 = oo PAM in a two-dimensional square lattice. Solid curves ( Vps and Vaf+fs ) represent the 
first-order phase transition and dotted curves (Vap) represent the second-order one. 




Fig. 2. V dependence of condensation energy AE, staggered magnetization of / electrons Mf and average /-electron density (nj). 
AE are variational energies of AFe (solid curves) and AFh (dotted curves) compared with that of PM. Open circles correspond to the 
transition points. 
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Fig. 3. Band dispersions of PM (y=1.50), AFh {V=1.00) and 
AFe (y=0.70) at Ep = -1.0 and n = 1.880. The Fermi energy 
is set to be 0. Arrows indicate the direction of the change of 
dispersion when V decreases (see text). 



It suggests the change of the topology of the Fermi sur- 
face and can be regarded as a kind of Fermi-surface re- 
construction. Although the condition of their calculation 
is somewhat different with ours (the chemical potential is 
fixed in ref. 19 instead of the electron density), their re- 
sult is consistent with the existence of the Fermi-surface 
reconstruction. 

Let us mention here the order of the transition. In 
our results, the Fermi-surface reconstruction is clearly a 
first-order transition. However, as mentioned before, it is 
to be noted that the second-order (or other continuous) 
transition cannot be completely ruled out since the result 
depends on the choice of the trial wave functions. 

Finally, we compare the present results with the Kondo 
lattice model (KLM). In the previous paper, we have 
shown that the Fermi-surface reconstruction occurs in 
the KLM. The result that the Fermi-surface reconstruc- 
tion occurs for [/ = cxd and sufficiently low value of E j in 
the PAM is consistent with the result of the KLM, since 
the KLM is an effective model of the PAM in the limit 
oiU ^ Qo and n/ ^ 1. Recent studies of the dynamical 
cluster approximation^*^ and the Gutzwiller approxima- 
tion^^ also support the existence of the Fermi-surface 
reconstruction in the KLM. In the Gutzwiller approxi- 
mation, the one-body part is more improved than ours 
but the obtained phase diagram is in good agreement 
with ours. This means that our variational wave func- 
tions capture the essence of the AF transition and the 
Fermi-surface reconstruction in spite of their rather sim- 
ple forms. 
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3.3 Momentum Distribution Function 

So far, we have discussed the Fermi surface obtained 
from the one-body part |$) without projection. However, 
it is in principle possible that the Fermi surface changes 
by the effect of the projection operator, Pq. Since the 
Fermi surface is determined from the discontinuities in 
the momentum distribution function, we calculate them 
for c and / electrons in this section. They are defined as, 

(8) 

(9) 

We calculate nc(fe) and n/(fc) with and without pro- 
jection, Pq, in order to study the role of the projec- 
tion operator. We set Ef = —1.0 and n = Uc + uf = 
612/324 = 1.889. Each k point is shifted along y direc- 
tion by (5 = 7r/9 because of the antiperiodic boundary 
condition. 

Figure 4 shows the result for PM a.tV = 1.50. In PM, 
the momentum distribution function of c and / electrons 
without projection are identical to u'l, (cq.(4)) and u| 
(eq.(5)). We denote them as n°(fc) and n"j{k), respec- 
tively. As shown by arrows in Fig. 4(b) and (c), n'^{k) 
and n° (fe) have discontinuities at the same position. Fig- 
ure 4(a) shows the Fermi surface determined from n°(fe) 
and n'j{k). We call it "original Fermi surface" in the fol- 
lowing. 

As shown in Fig. 4(b), the projection operator flattens 
n/(fc) in the whole Brillouin zone and the discontinu- 
ities are fairly reduced. Moreover, nf{k) has finite values 
even outside the original Fermi surface. They are known 
as a typical role of the on-site Coulomb interaction U. 
Judging from the fe-dependence of nf{k) near the orig- 
inal Fermi surface, it is probable that the position of 
the discontinuities does not change in the presence of 
Pq. Contrary to nf{k), ndk) is generally enhanced by 
projection since the excluded / electrons come into the 
c-clcctron orbital. As shown in Fig. 4(c), the discontinu- 
ities in ric(fc) are rather uncertain, but when we fit the 
data points by polynomials, we find that there are finite 
jumps at the original Fermi surface. Although it cannot 
be completely excluded that the discontinuities in nc{k) 
disappear, we expect that the original Fermi surface is 
preserved. If this is the case, we can say that our trial 
wave function satisfies the Luttinger's theorem: the vol- 
ume of the Fermi surface is determined by the electron 
density and not affected by the strength of electron cor- 
relation, unless the phase transition occurs. 

Next we discuss the AF state. As in the case of PM, 
n"(fc) and n"(fc) are explicitly calculated from the one- 
body part 1$). By diagonalizing the 4x4 matrix in 
eq.(7), we obtain four quasiparticles {otkai Pkailkai^ka) 
and band dispersions {E^,E'^, EI,EI) in the folded AF 
Brillouin zone. In the present condition, the Fermi energy 
crosses only the /3-band (see Figs. 3(b) and (c)). The a- 
band is completely filled while the 7- and the (5-band are 
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Fig. 4. Momentum distribution functions of PM at Ep = —1.0, 
n = 1.889 and V = 1.50. (a) Occupation of quasiparticles and 
Fermi surface in the first Brillouin zone. Solid (open) circles 
represent the occupied (unoccupied) fe-points. Each k point is 
shifted along y direction by 5 = 7r/9 because of the antiperiodic 
boundary condition. Calculated momentum distribution func- 
tions of (b) / electrons and (c) c electrons with and without 
C) projection. Arrows indicate the positions of discontinuities 
in n^(fe) and n^{k). The optimized variational parameters are 
(y,E/,Aj:,Ac) = (1.0899, -0.0460,0.0000, 0.0000) 



empty, i.e. 



i^, = 1 and 



nS - 



for all values of k. 



Q) = 



(10) 



Then n^{k) and n^^{k) are expressed as 

Ai{k) + 9{eF - E^)A5{k) 
= A2{k)+e{eF-E^)Ae{k) 
A3{k)+e{eF-E^)A^ik) 
U/(fe + Q) =Ai{k) + e{eF-E^)As{k) 

where 6{x) is the step function and ep denotes the Fermi 
energy. Here Ai{k) arc derived from the eigenvectors ob- 
tained in the diagonalization. We can see that the dis- 
continuities in n°(fe) and n^(fe) are determined by the 
occupation of the /3-band. 

Figure 5 shows the result for AFh at ^ = 1.00. 
The original Fermi surface determined from n°(fe) and 
n°(fc) is shown in Fig. 5(a). To make the comparison 
with the case of PM easier, we use the extended Bril- 
louin zone scheme in the following. The shaded areas in 
Figs. 5(b), (c) and (d) arc folded in the reduced Brillouin 
zone scheme. Compared with the case of PM, additional 




(71,71+8) 
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Fig. 5. Momentum distribution functions of AFh a-t Ep = —1.0, 
n = 1.889 and V = 1.00. (a) Occupation of quasiparticles in the 
/3-band and Fermi surface in the folded AF Brillouin zone. Solid 
(open) circles represent the occupied (unoccupied) fe-points. (b) 
First Brillouin zone described with the extended zone scheme. 
Calculated momentum distribution functions of (c) / electrons 
and (d) c electrons with and without (") projection. Arrows 
indicate the positions of discontinuities in n'j. (fc) and nj!(fe). 

The optimized variational parameters are (V, By, Ay, Ac) = 
(0.6368, -0.3236, 0.1193, 0.0000) 



discontinuities, whose magnitudes are Ar{k), appear in 
n°(fe) as shown in Fig. 5(c). It is originated from the ap- 
pearance of the AF order in / electrons. In the same way, 
additional discontinuities, whose magnitudes are A^lk), 
also appear in n°(fc), although they are very small. 

As shown in Fig. 5(c), nf{k) is flattened in the whole 
Brillouin zone by the effect of Pq as in the case of PM. 
However, discontinuities originated from nP^{k) seems to 
be preserved in nj{k). The discontinuities in nc{k) is 
rather uncertain because of their smallness (see the inset 
of Fig. 5(d)). However, by fitting the data points by poly- 
nomials, wc can find a finite jump at the original Fermi 
surface. In conclusion, we expect that the original Fermi 
surface is preserved through the projection also in AFh- 

Figure 6 shows the result for AFo at V = 0.90. Con- 
trary to the case of PM and AFh, the result is some- 
what different in AFe. The discontinuities exist near the 
AF Brillouin zone boundary both in n'^ik) and nPf{k) 
as shown by the arrows in Figs. 6(c) and (d). However, 
nf{k) is greatly flattened by projection and the discon- 
tinuities are greatly suppressed. If these discontinuities 
disappear, / electrons become completely localized and 
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Fig. 6. Momentum distribution functions of AFe at Ep = —1.0, 
n = 1.889 and V = 0.90. Notation is the same with that of Fig. 5. 
The optimized variational parameters are (y, Ay, Ac) = 
(0.6607, -0.0905, 0.4722, 0.0000) 



do not contribute to the Fermi surface. We cannot say 
whether the discontinuities in n/(fe) really disappear or 
not within this calculation. However, it is certain that 
the itinerancy of / electrons is greatly suppressed in AFe 
compared with the case of PM and AFh- On the other 
hand, the discontinuities in nc(fe) seems to exist at the 
original Fermi surface, as in the case of PM and AFh- 
In this sense, the volume and the shape of the Fermi 
surface is preserved through the projection also in AFe, 
whichever the discontinuities in nf{k) disappear or not. 
This result shows that even if the / electrons do not con- 
tribute to the Fermi surface, the c electrons "recognize" 
the existence of / electrons and preserve the original 
Fermi surface. It is possible that c electrons are redis- 
tributed by the effect of projection and form a "small" 
Fermi surface whose volume is dcitermined by the num- 
ber of c electrons alone. If it occurs, the discontinuities 
in nc(fc) exist only in the unshaded area and the discon- 
tinuities in the shaded area must disappear. But this is 
not the case. 

Finally we discuss how much the electrons are corre- 
lated and become heavy from the discontinuities in nc{k) 
and nf{k). We estimate the discontinuities at the Fermi 
surface from the difference between the values just in- 
side and just outside of the Fermi surface. As typical 
cases, we calculate them along two directions shown in 
the bottom panels in Fig. 7 and average their values. The 




Fig. 7. V dependences of the discontinuities of ric (fe) and nf{k) 
at the Fermi surface. They are estimated from the average of two 
directions shown above. 



upper panel of Fig. 7 shows the obtained Anc(fcF) and 
An/(A;F)- In PM and AFh, i^ricik-p) and An/(fcF) contin- 
uously decrease as V decreases. When the Fermi-surface 
reconstruction occurs (AFh AFo), a drastic change 
is observed. Anc(A:F) discontinuously increases since the 
position of ^f shifts near the AF Brillouin zone bound- 
ary, where the value of nc(fe) is rather large compared 
with the former position. It shows that the c electrons re- 
cover the weight at the Fermi energy. Note that even after 
the Fermi-surface reconstruction, c electrons are rcnor- 
malized compared with the case of AFs (Anc(fcF) = 1). 
On the other hand. An / (fcp) discontinuously decreases to 
a small value through the Fermi-surface reconstruction. 
It shows that the / electrons become almost localized 
and their itinerancy is fairly suppressed in AFe. From 
these results, we expect that the transition to AFe is a 
kind of orbital-selective Mott transition where only the / 
electrons become localized, as proposed in the CDMFT 
study. However, as mentioned before, we cannot say 
whether An/(A;F) becomes zero or not within this cal- 
culation. Moreover, the calculation for all directions is 
necessary for a precise discussion. It is left for an inter- 
esting future problem. 

4. Discussion 

We discuss the relation between our results and the 
experiments. If the Fermi-surface reconstruction occurs, 
changes in various physical quantities are expected. The 
abrupt change of dHvA frequencies across the AF tran- 
sition shows the appearance of a Fermi surface with dif- 
ferent topology, indicating the occurrence of a Fermi- 
surface reconstruction. The picture that the / electrons 
are localized in the AF phase and are itinerant in the 
paramagnetic phase, which is proposed from the com- 
parison of dHvA results and the band calculation, is also 
consistent with the scenario of a Fermi-surface recon- 
struction. 

On the other hand, it have been proposed that the AF 
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transition (x 0.8) and the abrupt change of dHvA fre- 
quencies {x ^ 0.4) are not coincident in CeRhi_a;CoxIn5. 
At X = 0.4, the magnetic structure changes from incom- 
mensurate {x < 0.4) to commensurate {x > 0.4) within 
the AF phase. ^ The fitting with band calculation sug- 
gests that the former AF phase corresponds to CeRhIn5 
type (4/-localized picture) and the latter corresponds to 
CeCoIn5 type (4/-itinerant picture). We expect that this 
transition is closely related to the Fermi-surface recon- 
struction within the AF phase (AF^ — >■ AFh), although 
the more realistic calculation is necessary for a precise 
discussion. 

The abrupt change of the Hall coefficient at the AF 
QCP in YbRh2Si2^ is also indicative of a Fermi-surface 
reconstruction. However, the band calculation shows that 
the Hall coefficient is determined by several contributions 
with opposite signs and is very sensitive to the energy 
level of / electrons. It was pointed out that small changes 
in the /-electron occupation are sufficient to reproduce 
the experimental result. Therefore, we cannot directly 
related the Fermi-surface reconstruction to the AF QCP 
in YbRh2Si2 at this stage. 

We also refer to the effective mass in the AF state. 
The dHvA experiment shows that the Fermi surface of 
CeRhIn5 is quite similar to that of non-4/ reference com- 
pound of LaRhlns.^^ It suggests that the contribution 
of the 4/ electrons to the Fermi surface is rather small 
in CeRhlns. However, the obtained cyclotron mass in 
CeRhIn5 is roughly by one order larger than that in 
LaRhIn5.^^ We consider that the existence of 4/ elec- 
trons does induce this mass enhancement, keeping the 
Fermi surface almost unchanged. As we have shown in 
the case of AFe, it is possible that / electrons do not 
contribute to the Fermi surface at all but the c electrons 
are renormalized and heavy reflecting the existence of the 
/ electrons. It can explain the enhancement of cyclotron 
mass in CeRhlns. 

5. Summary 

In summary, we have studied the ground state phase 
diagrams and some physical quantities of the PAM in a 
two-dimensional square lattice with VMC method. We 

have shown that 

• The Fermi-surfacc reconstruction which accompa- 
nies the localization of / electrons occurs in a wide 
region of parameters. It induces the topological 
change of the Fermi surface and stabilizes the lo- 
calized AF state. 

• The conventional itinerant AF transition also occurs 
separately from the Fermi-surface reconstruction. 

• Even in the localized AF state, the c electrons 
are renormalized compared with the non-interacting 
case, reflecting the existence of / electrons. 

These results will be a good starting point to clarify the 
relation between the AF transition and the change of 
the Fermi surface in heavy-fermion systems. For a more 



detailed discussion, calculations with realistic conditions 
will be necessary: three dimensional lattice structure, de- 
generacy of orbitals, magnetic ordering with incommen- 
surate wave vector, and so on. However, we think that the 
simple model discussed here capture the essence of the 
Fermi-surface reconstruction, which occurs in general as 
the evidence of strong localized character of / electrons. 
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